Teaching an Old Dog New Tricks: Improved Estimation of the Parameters of Stochastic Differential Equations by Numerical Solution of the Fokker-Planck Equation

Many stochastic differential equations (SDEs) do not have readily available closed-form expressions for their transitional probability density functions (PDFs). As a result, a large number of competing estimation approaches have been proposed in order to obtain maximum-likelihood estimates of their parameters. Arguably the most straightforward of these is one in which the required estimates of the transitional PDF are obtained by numerical solution of the Fokker-Planck(or forward-Kolmogorov) partial differential equation. Despite the fact that this method produces accurate estimates and is completely generic, it has not proved popular in the applied literature. Perhaps this is attributable to the fact that this approach requires repeated solution of a parabolic partial differential equation to obtain the transitional PDF and is therefore computationally quite expensive. In this paper, three avenues for improving the reliability and speed of this estimation method are introduced and explored in the context of estimating the parameters of the popular Cox-Ingersoll-Ross and Ornstein-Uhlenbeck models. The recommended algorithm that emerges from this investigation is seen to offer substantial gains in reliability and computational time.stochastic di®erential equations, maximum likelihood, ¯nite di®erence, ¯nite element, cumulative

Developing analytical distributions for temperature indices for the purposes of pricing temperature-based weather derivatives

Temperature-based weather derivatives are written on an index which is normally defined to be a nonlinear function of average daily temperatures. Recent empirical work has demonstrated the usefulness of simple time-series models of temperature for estimating the payoffs to these instruments. This paper develops analytical distributions of temperature indices on which temperature derivatives are written. If deviations of daily temperature from its expected value is modelled as an Ornstein-Uhlenbeck process with time-varying variance, then the distributions of the temperature index on which the derivative is written is the sum of truncated, correlated Gaussian deviates. The key result of this paper is to provide an analytical approximation to the distribution of this sum, thus allowing the accurate computation of payoffs without the need for any simulation. A data set comprising average daily temperature spanning over a hundred years for four Australian cities is used to demonstrate the efficacy of this approach for estimating the payoffs to temperature derivatives. It is demonstrated that expected payoffs computed directly from historical records is a particulary poor approach to the problem when there are trends in underlying average daily temperature. It is shown that the proposed analytical approach is superior to historical pricing.Weather Derivatives, Temperature Models, Cooling Degree Days, Maximum Likelihood Estimation, Distribution for Correlated Variables

Estimating the Payoffs of Temperature-based Weather Derivatives

Temperature-based weather derivatives are written on an index which is normally defined to be a nonlinear function of average daily temperatures. Recent empirical work has demonstrated the usefulness of simple time-series models of temperature for estimating the payoffs to these instruments. This paper argues that a more direct and parsimonious approach is to model the time-series behaviour of the index itself, provided a sufficiently rich supply of historical data is available. A data set comprising average daily temperature spanning over a hundred years for four Australian cities is assembled. The data is then used to compare the actual payoffs of temperature-based European call options with the expected payoffs computed from historical temperature records and two time-series approaches. It is concluded that expected payoffs computed directly from historical records perform poorly by comparison with the expected payoffs generated by means of competing time-series models. It is also found that modeling the relevant temperature index directly is superior to modeling average daily temperatures.Temperature, Weather Derivatives, Cooling Degree Days, Time-series Models.

The Devil is in the Detail: Hints for Practical Optimisation

Finding the minimum of an objective function, such as a least squares or negative log-likelihood function, with respect to the unknown model parameters is a problem often encountered in econometrics. Consequently, students of econometrics and applied econometricians are usually well-grounded in the broad differences between the numerical procedures employed to solve these problems. Often, however, relatively little time is given to understanding the practical subtleties of implementing these schemes when faced with illbehaved problems. This paper addresses some of the details involved in practical optimisation, such as dealing with constraints on the parameters, specifying starting values, termination criteria and analytical gradients, and illustrates some of the general ideas with several instructive examples.gradient algorithms, unconstrained optimisation, generalised method of moments.

Discrete time-series models when counts are unobservable

Count data in economics have traditionally been modeled by means of integer-valued autoregressive models. Consequently, the estimation of the parameters of these models and their asymptotic properties have been well documented in the literature. The models comprise a description of the survival of counts generally in terms of a binomial thinning process and an independent arrivals process usually specified in terms of a Poisson distribution. This paper extends the existing class of models to encompass situations in which counts are latent and all that is observed is the presence or absence of counts. This is a potentially important modification as many interesting economic phenomena may have a natural interpretation as a series of 'events' that are driven by an underlying count process which is unobserved. Arrivals of the latent counts are modeled either in terms of the Poisson distribution, where multiple counts may arrive in the sampling interval, or in terms of the Bernoulli distribution, where only one new arrival is allowed in the same sampling interval. The models with latent counts are then applied in two practical illustrations, namely, modeling volatility in financial markets as a function of unobservable 'news' and abnormal price spikes in electricity markets being driven by latent 'stress'.Integer-valued autoregression, Poisson distribution, Bernoulli distribution, latent factors, maximum likelihood estimation

Study of porous wall low density wind tunnel diffusers

Porous wall wind tunnel diffusers used with low density hypersonic nozzl

Point and interval estimation in two-stage adaptive designs with time to event data and biomarker-driven subpopulation selection

In personalized medicine, it is often desired to determine if all patients or only a subset of them benefit from a treatment. We consider estimation in two‐stage adaptive designs that in stage 1 recruit patients from the full population. In stage 2, patient recruitment is restricted to the part of the population, which, based on stage 1 data, benefits from the experimental treatment. Existing estimators, which adjust for using stage 1 data for selecting the part of the population from which stage 2 patients are recruited, as well as for the confirmatory analysis after stage 2, do not consider time to event patient outcomes. In this work, for time to event data, we have derived a new asymptotically unbiased estimator for the log hazard ratio and a new interval estimator with good coverage probabilities and probabilities that the upper bounds are below the true values. The estimators are appropriate for several selection rules that are based on a single or multiple biomarkers, which can be categorical or continuous